Theoretical study of the insulating oxides and nitrides: 
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An extensive theoretical study is performed for wide bandgap crystalline oxides and nitrides, 
namely, SiC>2, Ge02, AI2O3, SisN4, and Ge3N4. Their important polymorphs are considered which 
are for Si02i Q-quartz, a- and /3-cristobalite and stishovite, for GeC>2: a-quartz, and rutile, for 
AI2O3: a-phase, for SiaN4 and Ge3N4: a- and /3-phases. This work constitutes a comprehensive 
account of both electronic structure and the elastic properties of these important insulating ox- 
ides and nitrides obtained with high accuracy based on density functional theory within the local 
density approximation. Two different norm-conserving ab initio pseudopotentials have been tested 
which agree in all respects with the only exception arising for the elastic properties of rutile Ge02- 
The agreement with experimental values, when available, are seen to be highly satisfactory. The 
uniformity and the well convergence of this approach enables an unbiased assessment of important 
physical parameters within each material and among different insulating oxide and nitrides. The 
computed static electric susceptibilities are observed to display a strong correlation with their mass 
densities. There is a marked discrepancy between the considered oxides and nitrides with the latter 
having sudden increase of density of states away from the respective band edges. This is expected to 
give rise to excessive carrier scattering which can practically preclude bulk impact ionization process 
in Si 3 N 4 and Ge 3 N 4 . 

PACS numbers: 61.50.Ah, 62.20.Dc, 71.20.-b, 71.20.Ps, 77.22.Ch 



-a 

G 

O 

Q : 

(N 
> 



O 



1 

C 

o 
o 



X 



I. INTRODUCTION 



Insulating oxides and nitrides are indispensable ma- 
terials for diverse applications due to their superior 
mechanical, thermal, chemical and other outstanding 
high temperature properties. Furthermore, in the elec- 
tronic industry these wide band gap materials are be- 
ing considered for alternative gate oxides 1 - and in the 
field of integrated optics they provide low-loss dielectric 
waveguides- 2 -. Recently the subject of wide bandgap ox- 
ides and nitrides have gained interest within the con- 
text of nanocrystals which offer silicon-based technology 
for light emitting devices and semiconductor memories 3 . 
These nanocrystals are embedded in an insulating ma- 
trix which is usually chosen to be silica^&i. However, 
other wide bandgap materials are also employed such 
as germania 8 ^, silicon nitrid e 10 ' 11 ' 12 , and alumin a 13 ' 14 ' 15 . 
As a matter of fact, the effect of different host matrices 
is an active research topic in this field. 



Among these insulating oxides and nitrides technolog- 
ically most important ones are SiC>2, AI2O3, SiaN4. The 
activity around GeC>2 is steadily increasing. Another 
closely-related material, Ge3N4 has attracted far less at- 
tention up to now even though it has certain interesting 
properties 1 ^. The major obstacle has been the sample 
growth. However, a very recent study reported an in 
situ Ge3N4 growth on Ge, demonstrating high thermal 
stability and large band offsets with respect to the Ge 
system 1 ^. In this comprehensive work, we present the ab 
initio structural and electronic properties of all these ma- 
terials considering their common polymorphs; these are 
for Si02: a-quartz, a- and /3-cristobalite and stishovite 



phases, for Ge02: a-quartz, and rutile phases, for SisN4 
and Ge3N4i a- and /3-phases and for AI2O3: a-phase. 
For amorphous and inherently imperfect matrices, these 
perfect crystalline phases serve as important reference 
systems. Moreover, due to their distinct advantages, epi- 
taxial host lattices are preferred over the amorphous ones 
for specific applications. 

With an eye on these technological applications, we fo- 
cus on several physical properties of these lattices. The 
elastic constants play an important role on the strain pro- 
file of the embedded core semiconductor. Using Eshelby's 
continuum elastic consideration 18 the radial and tangen- 
tial stress fields of the nanocrystal can be determined 1 ^; 
these in turn, affect the optical properties-. The static 
and optical dielectric constants of these lattices intro- 
duce nontrivial local field effects that modify the absorp- 
tion spectra of an isolated nanocrystal when embedded 
inside one of these matrices^ -. Based on the simple ef- 
fective medium theory which has been tested by ab ini- 
tio calculations 21 -, one can assess which host lattice and 
nanocrystal combination would possess the desired op- 
tical properties. Because of the dielectric mismatch be- 
tween the nanocrystal core and the surrounding lattice, 
image charges will be produced 2 ^. These image charges 
should be taken into account in characterizing nanocrys- 
tal excitons22,. Another promising application is the visi- 
ble and near infrared electroluminescence from Si and Ge 
nanocrystals^. The electroluminescence is believed to be 
achieved by the recombination of the electron hole pairs 
injected to nanocrystals under high bias 3 . In this con- 
text the bulk state impact ionization process which can 
also give rise to electroluminescence is considered to be 
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detrimental leading to dielectric breakdown. For high- 
field carrier transport, the crucial physical quantity was 
identified to be the valence and conduction band density 
of states (DOS) for each of the crystalline polymorph 24 . 
Based on these technology-driven requirements we com- 
pute the elastic constants, band structures, dielectric per- 
mittivities and electronic DOS of these aforementioned 
crystal polymorphs. Our ab initio framework is based on 
the density functional theor y 25 ' 26 , using pseudopotentials 
and a plane wave basis 2 ^. With the exception of Ge 3 N 4 
which was far less studied, vast amount of theoretical 
work is already available spread throughout the literature 
based on a variety of technique8^' 29 i 30 t 31 ' 32 » 33 i 34 » 35 i 36 ^I. 
Our first-principles study here enables a uniform compar- 
ison of important physical parameters within each mate- 
rial and among different insulating oxides and nitrides. 

The plan of the paper is as follows: in Sec. II we provide 
details of our ab initio computations, Sec. Ill contains 
our first-principles results for the structural, electronic 
properties of the materials considered followed by our 
conclusions in Sec. IV. 



II. DETAILS OF AB INITIO COMPUTATIONS 

Structural and electronic properties of the polymorphs 
under consideration have been calculated within the den- 
sity functional theor y 25 ' 26 , using the plane wave basis 
pseudopotential method as implemented in the ABINIT 
code 2 ^. The results are obtained under the local density 
approximation (LDA) where for the exchange-correlation 
interactions we use the Tetcr Padc parameterization 38 , 
which reproduces Perdew-Zunger 39 - (which reproduces 
the quantum Monte Carlo electron gas data of Ceperley 
and Alder—). We tested the results under two differ- 
ent norm-conserving Troullier and Martins 41 " type pseu- 
dopotentials, which were generated by A. Khein and 
D.C. Allan (KA) and Fritz Haber Institute (FHI). For 
both pseudopotentials, the valence configurations of the 
constituent atoms were chosen as N(2s 2 p 3 ), 0(2s 2 p 4 ), 
Al(3s 2 3p 1 ), Si(3s 2 3p 2 ), and Ge(4s 2 4p 2 ). The number of 
angular momenta of the KA (FHI) pseudopotentials and 
the chosen local channel were respectively, for N: 1, p 
(3, d), for O: 1, p (3, d), for Al: 2, d (3, d), for Si: 2, 
d (3, d), and for Ge: 1, p (3, s). Our calculated val- 
ues for these two types of pseudopotentials were very 
similar, the only exceptional case being the elastic con- 
stants for rutile Ge02- Dielectric permitivity and the 
fourth-order tensor of elastic constants of each crystal 
are determined by starting from relaxed unit cell under 
the application of finite deformations within density func- 
tional perturbation theory 4 ^, as implemented in ABINIT 
and ANADDB extension of it. Another technical de- 
tail is related with the element and angular momentum- 
resolved partial density of states (PDOS). To get a rep- 
resentative PDOS behavior we need to specify the spher- 
ical regions situated around each relevant atomic site. 
The radii of these spheres are chosen to partition the 



bond length in proportion to the covalent radii of the 
constituent atoms. This resulted in the following radii: 
for the a-quartz Si02, rsi = 0.97 A, ro = 0.65 A, for 
the rutile Ge02, rg e = 1.16 A, ro = 0.69 A, for the a- 
A1 2 3 , r A i = 1.32 A, r = 0.56 A, and for the /?-Si 3 N 4 , 
rsi — 1.03 A, — 0.70 A. It should be pointed that 
even though such an approach presents a good relative 
weight of the elements and angular momentum channels, 
it inevitably underestimates the total DOS, especially for 
the conduction bands. Other details of the computations 
are deferred to the discussion of each crystal polymorph. 



III. FIRST-PRINCIPLES RESULTS 

First, we address the general organization and the un- 
derlying trends of our results. The lattice constants and 
other structural informations of all crystals are listed in 
Table fl] Table [IT] contains the bond lengths and bond 
angles of the optimized oxide polymorphs. These re- 
sults can be used to identify the representation of each 
polymorph within the amorphous oxides^ 2 -. The elas- 
tic constants and dielectric permittivity tensor of each 
crystal are tabulated in Table IIIII and Table IIV1 respec- 
tively. Very close agreement with the existing experi- 
mental data and previous calculations can be observed 
which gives us confidence about the accuracy and con- 
vergence of our work. Employing KA pseudopotentials, 
the band structure for the crystals are displayed along the 
high-symmetry lines in Figs. [TT3I5I6I together with their 
corresponding total DOS. Such an information is parti- 
clulary useful in the context of high-field carrier trans- 
port. These results are in good agreement with the pre- 
vious computation s 29 ' 32 ' 35 ' 36 . For all of the considered 
polymorphs the conduction band minima occur at the T 
point whereas the valence band maxima shift away from 
this point for some of the phases making them indirect 
band gap matrices (see Table IVf . However, the direct 
band gap values are only marginally above the indirect 
band gap values. These LDA band gaps are underesti- 
mated which is a renown artifact of LDA for semiconduc- 
tors and insulators^ 9 -. In this work we do not attempt any 
correction procedure to adjust the LDA band gap values. 

We present in Figs. I2I4I7I the element- and angular 
momentum-resolved PDOS. A common trend that can 
be observed in these various lattices is that their valence 
band maxima are dominated by the p states belonging to 
O atoms; in the case of SisN4 and GesN4 they are the N 
atoms. For the conduction band edges, both constituent 
elements have comparable contribution. This parallels 
the observation in amorphous Si02 where due to large 
electronegativity difference between Si and O, the bond- 
ing orbitals have a large weight on O atoms whereas the 
lowest conduction band states with antibonding charac- 
ter have a significant contribution from the Si atoms 60 . 

From another perspective, the band structures and the 
associated DOS reveal that there is a marked discrepancy 
between the valence and conduction band edges where for 
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TABLE I: Structural information on crystals. 



V_yl y a Lai 


vyi y aidi 
Structure 
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2.698 
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Tpt con fi 1 
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4 


2.372 


/3-cris. Si02 


Cubic 


a —7 403° 7 330° 7 160° 


Fd3m 


2 


1.966 


Stishovite Si02 


Tetragonal 


a =4.175° 4.145" 4.179 d 

C — Z.ODZ Z.04O Z.DDO 




2 


4.298 


a-quartz Ge02 


Hexagonal 


a =4.870° 4.861* 4.984 ; 

C — 0.0O4 O.OZU O.ODU 


P3 2 21 


3 


4.612 


Rutile Ge0 2 


Tetragonal 
Tetragonal 


a =4.283° 4.314" 4.4066° 
C=2.(8z 2.804 2.8619 


P42/mnm 


2 


6.655 


Q-AI2O3 


Rombohedral 


a =4.758° 4.762° 
c =12. ye 12.89b 


R3c 


2 


3.992 


a-Si 3 N 4 


Hexagonal 


a =7.732° 7.766' 
c =5.603° 5.615' 




4 


3.211 


/?-Si 3 N 4 


Hexagonal 


a =7.580° 7.585 J 
c =2.899° 2.895 J 


C'ih 


2 


3.229 


a-Ge 3 N 4 


Hexagonal 


a =7.985 a 
c =5.786 a 


r 4 


4 


5.691 


/3-Ge 3 N 4 


Hexagonal 


a =7.826 a 
c =2.993 a 




2 


5.727 



"This Work KA 
6 This Work FHI 
c Ref.|I3 
d Ref. 
fRei. 
SRef. 
e Ref. 
'Ref. 
JRef. 



TABLE II: Bond lengths and bond 


angles (in 


degrees) of SiO 


2 and GeO 


2 polymorphs where 


x represents a Si or a 


Ge atom. 


Crystal 




x-o (A) 


x-0 (A) 


O-x-0 


O-x-O 


O-x-0 


O-x-O 


x-O-x 


x-O-x 


a-quartz Si02 


This Work 
Exp.° 


1.613 
1.605 


1.618 
1.614 


110.75 
110.50 


109.32 
109.20 


109.07 
109.00 


108.47 
108.80 


140.55 
143.7 




a-quartz Ge02 


This Work 


1.693 


1.699 


113.03 


110.62 


107.94 


106.16 


130.56 




a-cris. Si0 2 


This Work 
Exp. 6 


1.597 
1.603 


1.596 
1.603 


111.59 
111.40 


110.08 
110.00 


109.03 
109.00 


108.02 
108.20 


146.02 
146.5 




/3-cris. Si0 2 


This Work 
Exp.° 


1.603 
1.611 




109.47 
107.80 








180 
180.00 




Stishovite Si02 


This Work 
Exp. d 


1.804 
1.760 


1.758 
1.810 


98.47 


81.53 






130.76 
130.60 


98.47 


Rutile Ge0 2 


This Work 


1.848 


1.824 


99.34 


80.66 






99.34 


130.33 



"Ref. |H 
6 Ref. SI 
°Ref.[l3 
°Ref.[5l] 



the former there occurs a sharp increase of DOS just be- 
low the band edge. As the probabilities of most scattering 
processes are directly proportional to DOS^i, in the case 
of high-field carrier transport the electrons should en- 
counter far less scatterings and hence gain much higher 
energy from the field compared to holes. In this respect 



Si3N 4 and Ge3N 4 are further different from the others 
where for both conduction and valence bands the DOS 
dramatically increases (cf. Fig. so that the carriers 
should suffer from excessive scatterings which practically 
precludes the bulk impact ionization for this material. 

Another common trend can be investigated between 



4 



TABLE III: Elastic constants and bulk modulus for each crystal. 



Crystal 


(GPa) 


Cn 


C12 


C13 


C14 


C33 


C44 


Ce6 


B 


a-quartz Si02 


KA 


76.2 


11.9 


11.2 


-17.0 


101.7 


54.0 


32.1 


35 




FHI 


79.5 


9.73 


9.54 


-18.9 


101.7 


55.5 


34.9 


35 




Exp. a 


87.0 


7.00 


13.0 


-18.0 


107.0 


57.0 


40.0 


38 




Exp. 6 


87.0 


7.00 


19.0 


-18.0 


106.0 


58.0 




40 


a-Cris. Si0 2 


KA 


49.30 


5.26 


-11.41 




44.78 


74.15 


26.85 


12 


/3-Cris. Si0 2 


KA 


194.0 


135.0 








82.67 




155 




FHI 


196.1 


134.2 








85.40 




155 


Stishovite Si02 


KA 


447.7 


211.0 


203.0 




776.0 


252.0 


302.0 


306 




FHI 


448.8 


211.1 


191.0 




752.0 


256.5 


323.0 


302 




Exp. c 


453.0 


211.0 


203.0 




776.0 


252.0 


302.0 


308 


a-quartz Ge02 


KA 


66.7 


24.3 


23.1 


-3.00 


118.7 


41.3 


21.2 


41 




FHI 


63.8 


25.7 


26.2 


-0.81 


120.2 


35.3 


19.1 


42 




Exp. a 


66.4 


21.3 


32.0 


-2.20 


118.0 


36.8 


22.5 


42 




Exp. 6 


64.0 


22.0 


32.0 


-2.00 


118.0 


37.0 


21.0 


42 


Rutile Ge0 2 


KA 


405.9 


235.3 


189.2 




672.4 


206.0 


314.4 


292 




FHI 


349.2 


197.2 


185.1 




617.5 


171.8 


274.8 


258 




Exp. c 


337.2 


188.2 


187.4 




599.4 


161.5 


258.4 


251 


Q-AI2O3 


KA 


493.0 


164.1 


130.1 




485.8 


155.5 


164.4 


258 




Exp/ 


497.0 


164.0 


111.0 




498.0 


147.0 




251 


/3-Si 3 N 4 


KA 


421.8 


197.8 


116.6 




550.7 


100.2 


112.0 


250 




Exp. 9 


433.0 


195.0 


127.0 




574.0 


108.0 


119.0 


259 




Exp.' 1 


439.2 


181.8 


149.9 




557.0 


114.4 


135.9 


265 


/3-Ge 3 N 4 


KA 


364.3 


184.9 


111.7 




486.3 


80.4 


89.7 


225 



a Ref.[H 
'Ref.m 
c Ref. 
d Ref. 
e Ref. 
•fRef. 
»Ref. 
ft Ref. 



TABLE IV: Dielectric permittivity tensor. 



Crystal 


c ix yy 


e u 


OO OO 

'zx yy 


OO 

c zz 


Q-quartz Si02 


4.643 


4.847 


2.514 


2.545 


a-cris. Si02 


4.140 


3.938 


2.274 


2.264 


/9-cris. Si02 


3.770 


3.770 


2.078 


2.078 


Stishovite Si02 


10.877 


8.645 


3.341 


3.510 


a-quartz Ge02 


5.424 


5.608 


2.864 


2.947 


Rutile Ge0 2 


10.876 


8.747 


3.679 


3.945 


Q-AI2O3 


10.372 


10.372 


3.188 


3.188 


0-Si 3 N 4 


8.053 


8.053 


4.211 


4.294 


0-Ge 3 N 4 


8.702 


8.643 


4.558 


4.667 



the density of each polymorph and the corresponding 
static permittivity, e s . Such a correlation was put for- 
ward by Xu and Ching among the SiC>2 polymorphs 29 . 
We extend this comparison to all structures considered in 
this work and rather use % e = e s — 1 which corresponds 
to electric susceptibility. It can be observed from Fig. [8] 
that the trend established by Si02 polymorphs is also 
followed by /3-Si3N 4 and a-A^O^. On the other hand, 
Ge-containing structures while possessing a similar trend 
among themselves, display a significant shift due to much 
higher mass of the this atom. This dependence on the 



TABLE V: Indirect (E g ) and direct (E g (T)) LDA Band Gaps 
for each crystal. 



Crystal 


VB Max. 


CB Min. 


E g (eV) E g (T) (eV) 


a-quartz Si02 


K 


r 


5.785 


6.073 


a-cris. Si02 


L 


r 


5.525 


5.525 


/9-cris. Si02 


L 


r 


5.317 


5.317 


Stishovite Si02 


L 


r 


5.606 


5.606 


a-quartz Ge02 


K 


r 


4.335 


4.434 


Rutile Ge0 2 


L 


r 


3.126 


3.126 


a-Al 2 3 


L 


r 


6.242 


6.242 


a-Si 3 N 4 


M 


r 


4.559 


4.621 


/3-Si 3 N 4 


A-F 


r 


4.146 


4.365 


a-Ge 3 N 4 


M 


r 


3.575 


3.632 


,3-Ge3N4 


A-r 


r 


3.447 


3.530 



atomic mass needs to be removed by finding a more suit- 
able physical quantity. We should mention that such a 
correlation does not exist between the volume per prim- 
itive cell of each phase and the static permittivity. After 
these general comments, now we concentrate on the re- 
sults of each lattice individually. 
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(a) a - Quartz SiO, 




T K H A T M L A I) 15 

DOS 



(c) p - QrisWbalifc SiO, 




L r k xw r o is 

DOS 



FIG. 1: LDA band structure and total DOS (electrons/eV cell) 
Si0 2 , and (d) stishovite Si02- 

A. Si0 2 

The a-quartz Si02 is one of the most studied poly- 
morphs as it is the stable phase at the ambient pres- 
sure and temperatur e 30 ' 34 , furthermore its short-range 
order is essentially the same as the amorphous SiC>2 60 . 
a-quartz SiC>2 has a hexagonal unit cell containing three 
SiC>2 molecules. A plane-wave basis set with an en- 
ergy cutoff of 60 Ha was used to expand the electronic 
wave functions at the special fc-point mesh generated by 
10x10x8 Monkhorst-Pack scheme^ 2 -. The band structure 
of a-quartz SiC>2 has been calculated by many authors 
(see, for instanc o 28 ' 29 ). Our calculated band structure 
and total DOS shown in Fig. QJa) are in agreement with 
the published studies^ 9 -. The indirect LDA band gap for 
this crystal is 5.785 eV from the valence band maximum 
at K to the conduction band minimum at T. The direct 
LDA band gap at T is slightly larger than the indirect 
LDA band gap as seen in Table fVl Calculated values of 



(b) ct- Cristobal ito SiO, 




rZRXTMXO 15 

DOS 



Id) Stastevite SiO, 




rZRXTMXO 10 

DOS 



of (a) Q-cristobalite S1O2, (b) a-quartz Si02, (c) /3-cristobalite 



the elastic constants and bulk modulus listed in Table [TtTI 
are in good agreement with the experiments. Apart from 
C12, the elastic constants are within 10% of the experi- 
mental values. The discrepancy in C12 can be explained 
by the fact that C12 is very soft and this type of deviation 
also exists among experiments which is also the case for 
C14. 

a-cristobalite Si02 has a tetragonal unit cell contain- 
ing four Si02 molecules. In the course of calculations 
an absolute energy convergence of 10~ 4 Ha was obtained 
by setting a high plane wave energy cutoff as 60 Ha and 
10x10x8 fc-point sampling. Figure [Ijb) shows the band 
structure of a-cristobalite Si02 with the 5.525 eV direct 
band gap at T. The bulk modulus of 12 GPa is the small- 
est among all the host lattice polymorphs considered in 
this work. 

Regarding /3-cristobalite, its actual structure is some- 
what controversial, as a number of different symme- 
tries have been proposed corresponding to space groups 
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(a) 



(h) 





Tend 


■ 111 ik 


Si 







-20 -15 -10 -5 5 10 15 20 
Energy (eV) 



-10 -5 5 10 15 20 

Energy (eV) 



FIG. 2: DOS of a-quartz S1O2 (a) Element-resolved; total, 
PDOS of Si, PDOS of O. (b) Angular momentum-resolved; 
Si s electrons, Si p electrons, Si d electrons (not visible at the 
same scale), O s electrons, O p electrons. 



Fd3m, I42d, and P2i32i. Recently, incorporating the 
quasiparticle corrections the tetragonal IA2d phase was 
identified to be energetically most stable^ 3 -. However, we 
work with the structure having the space group of Fdim 
that was originally proposed by Wyckofl* 4 and which is 
widely studied primarily due to its simplicit y 28 ^ 30 . This 
phase has a cubic conventional cell with two molecules. 
We used 60 Ha plane wave energy cutoff and 10x10x10 
fc-point sampling. Figure [TJc) shows the band structure 
of /3-cristobalite SiC>2 with the 5.317 eV direct band gap 
at r. Unlike their band structures, total DOS of a- and 
and /3-cristobalite SiC>2 are very similar (cf. Fig. HJc)). 
This similarity can be explained by the fact that their 
local structure are very close. On the other hand there is 
a considerable difference between the DOS spectra of the 
a-quartz Si02 and the /3-cristobalite Si02- In Table HTT1 
we present elastic constants of the /3-cristobalite Si02 cal- 
culated by two types of pseudopotentials, FHI and KA. 
There is no considerable difference between them. Di- 
electric constants of /3-cristobalite Si02 are the smallest 
among the five polymorphs of Si02 studied here (see Ta- 
bleHVl). 



Stishovite is a dense polymorph of Si02 with octahe- 
drally coordinated silicon, unlike the previous phases 3 ^. 
It has a tetragonal cell with two molecules. Calcula- 
tions were done by using 60 Ha plane wave energy cutoff 
and 8x8x10 /fc-point sampling. The band structure of 
stishovite with a wide single valence band is markedly dif- 
ferent from that of the previous three crystalline phases of 
Si02 having two narrow upper valence bands. The cause 
of this increased valence bandwidth is the lack of separa- 
tion between bonding and nonbonding states 3 ^. Hence, 
the total DOS for stishovite shows no gap at the mid- 
dle of the valence band (see Fig.QJd)). Our calculations 
yield a direct LDA band gap of 5.606 eV at T. As seen in 
Table IIII1 the differences between our computed elastic 
constants and the experimental values are less than 3%; 
this is an excellent agreement for LDA. Its bulk modu- 



lus is the largest among all the host lattice polymorphs 
considered in this work. Moreover, dielectric constants 
of stishovite is the largest of the five polymorphs of Si02 
considered in this work (see Table HV|) . 



B. GeO z 

For a-quartz Ge02 we used the same energy cutoff and 
fc-point sampling as with a-quartz Si02 which yields ex- 
cellent convergence. The band structure of the a-quartz 
Ge02 is displayed in Fig. EJa). The similarity of the 
band structures of the a-quartz Ge02 and the a-quartz 
Si02 is not surprising as they are isostructural. Similarly 
their total DOS resemble each other (cf. Fig.[3Ja)). The 
indirect LDA band gap for this phase is 4.335 eV from 
the valence band maximum at K to the conduction band 
minimum at T . The direct band gap at T is slightly dif- 
ferent from indirect band gap as seen in Table [V] This 
gap is smaller than that of the a-quartz Si02- The per- 
fect agreement between calculated elastic constants of 
the a-quartz Ge02 and experimental value s 45 ' 54 can be 
observed in Table HTT1 

The rutile structure of Ge02 , also known as argutite^ 5 - 
is isostructural with the stishovite phase of Si02- The 
same energy cutoff and fc-point sampling values as for 
stishovite yield excellent convergence. The direct LDA 
band gap at T for rutile-Ge02 is less than that of 
stishovite with a value of 3.126 eV. The two upper va- 
lence bands are merged in the total DOS (see Fig. O^b)) 
as in the case of stishovite. The increased valence band- 
width in the band structure can be explained by the same 
reason as in the case of stishovite. The results of the 
elastic constants calculated with KA type pseudopoten- 
tial shown in Tabic IIIII deviate substantially from the 
experiment whereas the agreement with the FHI pseu- 
dopotentials is highly satisfactory. The similarity of the 
dielectric constants of rutile Ge02 and stishovite can be 
observed in Table HVl 



C. AI2O3 

AI2O3 is regarded as a technologically important oxide 
due to its high dielectric constant and being reasonably 
a good glass former after SiO^. The a-A^Os (sapphire) 
has the rhombohcdral cell with two molecules. Computa- 
tions about AI2O3 were done by using 60 Ha plane wave 
energy cutoff and a total of 60 fc-points within the Bril- 
louin zone. Fig. O shows the computed band structure 
and total DOS of the a-A^Os. These are in excellent 
agreement with the previous calculatio n 31 ! 33 . For AI2O3, 
minimum of the conduction band is at T and maximum 
of the valence band is at a point along T — X close to the 
r point. The corresponding LDA band gap is 6.242 eV. 
Because of the very small difference between the direct 
and indirect band gaps, AI2O3 is considered as a direct 
band gap insulator. Measured band gap of this crystal 
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FIG. 3: LDA band structure and total DOS of (a) a-quartz GeC>2, (c) rutile GeCh- 



is 8.7 eV. However the precise value of the gap of AI2O3 
is still elusive because of the existence of an excitonic 
peak near the absorbtions edge^. As seen in Table IIII) 
computed values of the elastic constant and bulk modu- 
lus of AI2O3 are in excellent agreement with the exper- 
iments. As a furher remark, the a-Al2 0,3 unit cell can 
be described as hexagonal or rhombohcdral depending 
on the crystallographical definition of the space group 
R3C. During our first-principles calculations it has been 
defined as rhombohedral in which case C14 vanishes. Al- 
though the sign of C14 is experimentally determined to be 
negative for the hexagonal- AI2O3, previous calculations 
reported a positive valued. To check this disagreement 
we have calculated the elastic constant of the hexagonal- 
AI2O3 and found it to be around -3.0. 
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FIG. 5: LDA band structure of and total DOS of a-Al 2 3 . 
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FIG. 4: DOS of rutile Ge0 2 (a) Element-resolved; total PDOS 
of Ge, PDOS of O. (b) Angular momentum-resolved; Ge s 
electrons, Ge p electrons, Ge d electrons, O s electrons, O p 
electrons. 



D. Si 3 N 4 and Ge 3 N 4 

The research on silicon nitride has largely been driven 
by its use in microelectronics technology to utilize it as an 
effective insulating material and also as diffusion mask for 
impurities. Recently it started to attract attention both 
as a host embedding material for nano cry st al o 1 ' 1 1 1 1 2 and 
also for optical waveguide applications 2 . The a- and j3- 
SisN4 have hexagonal conventional cells with four and 
two molecules, respectively. We used 60 Ha plane wave 
energy cutoff and 6x6x8 fc-point sampling. The com- 
puted band structures of these two phases shown in 
Figs. [B] (a) and (b) are identical to those reported by 
Xu and Ching 32 . The top of the valence band for (3- 
Si3N4 is along the T-A direction, and for a-Si3N 4 it is 
at the M point. The bottom of the conduction band 
for two phases are at the T point. The direct and indi- 
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FIG. 6: LDA band structure and total DOS of (a) a-Si 3 N 4 , (b) /3-Si 3 N 4 , (c) a-Ge 3 N 4 and (d) /3-Ge 3 N 4 . 
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FIG. 7: Element-resolved DOS of (a) /3-Si 3 N 4 ; total, PDOS 
of Si, PDOS of N, (b) /3-Ge 3 N 4 ; total, PDOS of Ge, PDOS of 
N. 



rect LDA band gaps of these two phases are respectively, 



4.559 eV, 4.621 eV for a-Si 3 N 4 and 4.146 eV, 4.365 eV for 
the /?-Si3N 4 . The general band structure of two phases 
are very similar, except that the a-Si3N 4 has twice as 
many bands because the unit cell is twice as large. The 
total DOS of these two phases shown in Figs. [6fa) and 
(b) are only marginally different. Calculated values of 
the elastic constants and bulk modulus of /3-Si3N 4 listed 
in Table IIIII are in excellent agreement with the quoted 
experiments. Those for the a-Si3N 4 which is thermody- 
namically less stable with respect to /3-phase^ were left 
out due to excessive memory requirements for the desired 
accuracy. 

Ge3N 4 is the least studied material among the ox- 
ides and nitrides considered in this work. Recently 
its high-pressure 7-phase has attracted some theoretical 
interest. However, the available Ge3N 4 samples contain 
a mixture of a and /3-phases as in the case of Si3N 4 and 
these are the polymorphs that we discuss in this work. 
The band structures of both of these phases of Ge3N 4 (cf. 
Fig. are very similar to those of S13N4. Regarding the 
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elastic constants of /3-Ge3N4, our theoretical results listed 
in Table IIIII await experimental verification. In terms of 
density, the (3 phases of SiaN4 and Ge3N4 fill the gap be- 
tween the a-quartz and stishovite/rutile phases of their 
oxides. As can be observed from Fig. [8] their electric sus- 
ceptibility versus density behavior strengthens the corre- 
lation established by the remaining polymorphs. Finally 
it should be pointed that /3-Ge3N4 has the largest high- 
frequency dielectric constant (eoo) among all the materi- 
als considered in this work. 



VsiO, 

8 - 

§ - 
4 - 

• ac-SiO, 

Pc-SrO, 



Density (gr/cm ) 



technologically-important insulating oxides and nitrides: 
Si0 2 , Ge0 2 , A1 2 3 , Si 3 N 4 , and Ge 3 N 4 . The structural 
parameters, elastic constants, static and optical dielec- 
tric constants are obtained in close agreement with the 
available results. The computed dielectric constants are 
observed to display a strong correlation with their mass 
densities. For all of the considered polymorphs the con- 
duction band minima occur at the T point whereas the 
valence band maxima shift away from this point for some 
of the phases making them indirect band gap matrices. 
However, the direct band gap values are only marginally 
above the indirect band gap values. The investigation 
of band structure and DOS data reveal that the holes in 
all polymorphs considered and the electrons for the case 
of SisN4 and Ge3N4 should suffer excessive scatterings 
under high applied field which will preclude bulk impact 
ionization for these carrier types and polymorphs. This 
can be especially important for applications vulnerable 
to dielectric breakdown. 
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